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Abstract 

By means of variational methods, in this paper, we estabhsh sharp existence 
results for solutions of the master equations governing 'fractional multiple vortices.' 
In the doubly periodic situation, the conditions for existence are both necessary 
and sufficient and give the upper bounds for the vortex numbers in terms of the size 
of the periodic cell domain. In the planar situation, there is no restriction on the 
vortex numbers. In both situations, the solutions are uniquely determined by the 
prescribed locations and the local winding numbers of the vortices. 
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1 Introduction 



Recall that the classical Abrikosov-Nielsen-Olesen vortices [1, 19, 29] are static solutions 
of the Abehan Higgs model governed by the Hamiltonian 

n^\F' + \\Dcf>\' + ^{\cl>\^-l)\ (1-1) 

where F = 91^2 — <92^i is the magnetic field induced from the real- valued gauge potential 
field A = (741,742), 4> is the complex- valued scalar Higgs field, D(() = V0 — iA0 ~ 
D2(f)) denotes the gauge-covariant derivative, A > is a coupling parameter, which 
satisfy the Euler-Lagrange equations of (1.1) given by 

- 2iA • V0 - i(V • A)0 = |A|20+^(|0|2-l)0, (1.2) 
AA-V(V-A) = -{^D(f)- (f)D^), (1.3) 

also known as the Ginzburg-Landau equations [16]. At the critical coupling A = 1 sep- 
arating type-1 and type-11 superconductivity, the equations (1.2) and (1.3) over are 
equivalent [43] to the self-dual or BPS (after the pioneering explorations of Bogomol'nyi 
[5] and Prasad-Sommerfield [30]) vortex equations 

Di(j)±iD2(j) = 0, (1.4) 
F± 1(10^-1) = 0. (1.5) 

As a consequence, the solutions are all characterized precisely by the zero distribution 
of the Higgs field [19, 42, 43] where superconductivity is absent. More precisely, the 
zeros of 0, say Pi, ■ ■ ■ ,PAr, which may be arbitrarily prescribed, are the centers of magnetic 
vortices. Solutions of such type are known as multiple vortices, and analytically (cf. [19]) 
they are governed by the master equation 

N 

Ai^ = e"-l + 47r^(5p,(x). (1.6) 

Although the equation (1.6) seems to be rather simple, it is not integrable [31] and its 
understanding requires techniques [19, 41, 45, 46] based on functional analysis. Besides, 
an equation of similar structure appears in conformal geometry [2, 8, 9, 10, 11] which has 
been an actively pursued subject in geometric analysis. 

Physically, the solutions of (1.6) give rise to thinly formed vortex tubes in a super- 
conductor, a phenomenon due to the Meissner effect, which is essential to a mechanism 
called monopole confinement [17, 34]. Those facts provide some of the earliest clues 
[25, 26, 28, 36, 38] to the great puzzle of quark confinement. Indeed, for quark confine- 
ment, Seiberg and Witten [32] showed how vortex tubes are generated in non-Abelian 
gauge field theory models through a dual Meissner effect and how quark confinement is 
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achieved through (color-charged) non-Abehan monopole condensation [17, 34]. Due to 
the interest and importance of non-Abehan vortices in monopole and quark confinement 
problems, a burst of studies have been carried out and a wide range of nonlinear vortex 
equations of rich features have been derived by theoretical physicists. Here we mention 
the seminal papers of Marshakov and Yung [27], Hanany and Tong [18], and Auzzi, Bolog- 
nesi, Evslin, Konishi, and Yung [3], and refer to the survey articles [14, 20, 21, 33, 44] 
and references therein for various relevant results in this active area. Also we indicate in 
[22, 23, 24] a series of existence and uniqueness results concerning the associated partial 
differential equation problems. 

In the present paper, we aim at establishing the existence and uniqueness theorems for 
solutions of the master equations governing 'fractional vortices' in the recent work of Eto, 
Fujimori, Gudnason, Konishi, Nagashima, Nitta, Ohashi, and Vinci, [13]. These equations 
may be regarded as extensions of the classical master equation (1.6) governing the Abelian 
Higgs vortices. We shall show that, unhke the equations arising in the Weinberg-Salam 
electroweak theory [4, 12, 35, 40] and the Chern-Simons theory [6, 39, 41, 46], the elegant 
structures of these master equations allow us to achieve a complete understanding of 
their solutions in terms of explicitly stated conditions. In the next section, we describe 
the master equations obtained in [13], to be studied here. We state our main existence and 
uniqueness results in §2 and in the subsequent sections, §3 - §5, we give the proofs. Our 
methods here are centered around a direct minimization procedure, an idea initiated in 
the recent work [22]. We note that, while direct minimization is convenient to implement 
for problems over a doubly periodic domain (due to compactness), it is not an easy task 
to carry on for problems over the full plane (due to loss of compactness). In §5, we 
elaborate on such a problem and overcome the difficulty by establishing some coercive 
estimates for the action functionals. We shall see that the positive and negative parts of 
a trial function need to be treated separately and the coercive lower bounds for them are 
rather different. We anticipate that this new method will be useful in dealing with other 
problems of similar difficulties. In §7 and §8, we study and describe various degenerate 
situations. 



2 The master equations and existence theorems 

The first master equation obtained in [13] is a scalar equation of the form 

N 

Au = A(m|A|'e"^" + V" - + 47r^ 5^^, (2.1) 

s=l 

where ^ > 0, {A, B) G is a constant vector which lies in the elliptical vacuum manifold 
(denoted by M) defined by the equation 

m\Af + n\B\^ (2.2) 

and the set of points {pi, ■ ■ ■ ,Pn}, where the p's are not necessarily distinct, describes 
the locations of N vortices as in the Abehan Higgs case, (1.6). There are two problems of 
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interest: the solutions of (2.1) defined over a doubly (gauge) periodic domain fl, originally 
formulated in [37] and elaborated in concrete situations in [6, 35, 45], and those over the 
full plane M^. In the latter situation, the vacuum manifold (2.2) leads us to the following 
natural boundary condition for u: 

u ^ as |x| — >■ oo. (2.3) 

Our existence theorem for (2.1) may be stated as follows. 

Theorem 2.1. Consider the multiple vortex equation (2.1) with arbitrarily prescribed 
points pi, ■■■ ,pn, and assume that (2.2) holds subject to the non- degeneracy condition 
m,n > 0, and ^ > 0. Then 

(i) a solution exists over a doubly periodic domain Q if and only if 

AttN 

^<m\; (2.4) 

(a) a solution over the full plane subject to the boundary condition (2.3) always 
exists and tends to zero exponentially fast; 

(Hi) in both cases, the solution u is unique and satisfies 

u<0, , . 

m|ylp||e"™ - l||i + n|5p||e"" - = ^N, ^ ' 

where || • ||i denotes the L^-norm over Q or respectively. 

The second and third master equations obtained in [13] are systems of two nonlinear 
equations of the form 

A^/i = Ai(m|Ape™"^ + |5|V^+"2 + |C|Vi-"2 -Ci) +47r (2.6) 
Au2 = A2(|5|Vi+'^^ - |C|V^-"^ -e2) + 47r52^P2.., (2.7) 

s=l 

where m = 1 or 2, respectively, Ai, A2, ^i, > are all fixed parameters, and {A, B, C) e 
satisfies 

m|A|^+ + |Cp = 6, |5p-|C|2 = 6. (2.8) 

The vanishing of either one of the parameters B and C makes the system of equations 
(2.6)-(2.7) degenerate into the Abelian model (1.6), as will be discussed in §7 and §8. 

Thus, to capture the true new features of (2.6)-(2.7), we shall discuss it first under 
the non-degeneracy conditions 

By^O and C 0. (2.9) 
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Notice that, for C 7^ 0, we subject the sets of vortex points to the condition 

{P2,l, ■ ■ ■ ,P2,N2} C {pi,l, • • • ,Pl,iVi}, (2.10) 

and impose the natural restriction 

N2 < Ni, (2.11) 

which ensure the regularity of the quantity e"^~"^. 

From the vacuum manifold (2.8), we see that for a solution (1*1,1*2) of (2.6)-(2.7) over 
M^, the natural boundary condition is expressed as follows: 

1*1,1*2^0 as ^ 00. (2.12) 

On the other hand, the vacuum manifold (2.8) plays no role as far as the solvability 
of the system of equations (2.6)-(2.7) is concerned under periodic boundary conditions 
with a given cell domain Q. In this case, we identify necessary and sufficient conditions 
for the existence of a periodic solution that impose upper bounds on the vortex numbers 
A^i and N2 (see (2.13) and (2.14) below), a feature common to other self-dual periodic 
vortex problems (cf. [41, 46]). Actually, the vacuum constant enters to ensure suitable 
pointwisc estimates on the solution, consistently with their physical interpretation. 

More precisely, we prove 

Theorem 2.2. Consider the master equations (2.6) and (2.7) for any integer m > 
subject to the non- degeneracy condition (2.9). Then 

(i) a solution over a doubly periodic domain Q exists if and only if the conditions 

47r(^ + ^) < (^1 + 6)1^^1, (2.13) 
4^(^-^) < (Ci- 6)1^^1, (2.14) 

hold simultaneously; 

(a) a solution over the full subject to the boundary condition (2.12) always exists 
provided (2.8) holds; 

(Hi) in both cases above, when a solution exists, it must be unique; 
(iv) furthermore, if we assume (2.8) together with the condition 

A2>Ai, 6>0, and m|A|2 < min |l, £| - 1^ |C|^ (2.15) 

then 

Ml±M2<0, Ml < 0, (2.16) 
and the following "quantization " property hold: 

m|AH|l-e'""i||i + 2|Sn|l-e"i+"2||i = 47r f ^ + (2.17) 

\ Ai A2 / 

m|^n|l-e™"i||i + 2|C|2||l-e"i-"2||i = 47r ( ^ - ^ ) , (2.18) 

\ Ai A2 / 

where \\ ||i denotes the L}-norm, over or'M?, accordingly. 
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Remark 2.3. The conditions (2.11) and (2.18) almost justify the necessity of the condi- 
tion A2 > Ai. 

For convenience, degenerate cases where B or C will be allowed to vanish will be 
treated separately in latter sections. 



3 Proof of Theorem 2.1 

Wc first consider (2.1) over a doubly periodic domain Vt. Let uq be the unique periodic 
solution over of the problem: 

/^Modx = 0. 



Thus, with u — uq-\-v, the equation (2.1) subject to periodic boundary condition becomes: 

(3.2) 



Aw = A(m|A|2e'""o+™'' + n|5| - + f^T' ^ ^ ^' 
V is doubly periodic over dO.. 

Integrating (3.2) over O, we have 

/" AttN 

/ (m|A|2e'""°+'"'' + n|5|2e""°+"'') dx = — = ?7 > 0, (3.3) 

JQ. A 

which gives us the condition (2.4) as a necessary condition for a solution to exist. We 
now show that (2.4) is sufficient to ensure the existence of a solution. 

It is clear that (3.2) is the Euler-Lagrange equation of the variational functional 

I[V) = J |^|Vi;|' + A (|^|2e™-0+m^ + |^|2gn«o+n.) _ ^^^^ _ ^| (3 4) 

over the Sobolev space 

H{fl) — {v & W^^^{M.'^) I V is doubly periodic with periodic cell domain fl}, (3.5) 

equipped with the standard scalar product and corresponding norm. To proceed further, 
we decompose H{Q) = R © H{Q) with 

f = L+f. /eK, /e Hip) satisfying [ f dx ^ 0. (3.6) 
By means of the Jensen inequality, we have: 

> ^ / iVi-I'dx + Aini (|A|'e"*^+ |S|V^) - A77T;. (3.7) 
2 Jn 



Since for a > 0, & > 0, a + 6 > 0, c > 0, m > 0, n > 0, and am + bn > 0, the function 

a{t) = ae"** + 6e"* - ct, (3.8) 

satisfies: a(t) — )■ oo as t — )■ ±oo, we see from (3.7) that /(■) is bounded from below and, 
by virtue of the Poincare inequahty, it is also coercive in H{Q). By the weak lower semi- 
continuity of / and its strict convexity, we conclude that / attains the minimum value at 
a point that is its unique critical point. So part (i) of Theorem 2.1 is established. 

Now consider (2.1) over the full plane. In view of the condition (2.2), it may be 
rewritten as follows: 

N 

Au = Am|A|2(e'"" - 1) + An|B|2(e"" - 1) + ^ttJ^^Ps- (3-9) 



We introduce the background function: 

N 



/I |2 \ 



satisfying: 

N 1 

Auo^AnJ2Sp.-9o, go^^l^Yl (^+|^_p^|2)2 - (3-11) 
Thus, setting u — Uq + v, we recast (3.9) into the equation 

Av = Am|A|2(e'""°+'"" - 1) + An|5|2(e'^'^°+™ - 1) + ^o, xe M^ (3.12) 
corresponding to the Euler-Lagrange equation of the functional 

I{v) = l-\\Vv\\l + X\Af I {e'""°(e"*^ -l-mv) + (e™"° - l)mv} dx 

+A|5|2 [ {e"""(e"'' -l-nv) + (e""° - l)nv} dx + [ govdx, (3.13) 

well defined and of class over the Sobolev space VF^'^(M^). Indeed, the above functional 
is similar to that derived for the equation (1.6) in [19, 42]. Hence, using the method in 
[19, 42], the existence and uniqueness of a critical point of / in the space W^''^{M.'^) can 
be established. Similarly, (2.5) follows for example by arguing as in Proposition 3.2.4 of 
[41]. Here, we omit the details. 

We note that the existence of a critical point of the functional (3.13) may also be 
established by our direct minimization method to be presented in §5. 
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4 Proof of Theorem 2.2: doubly periodic case 

We now consider the equations (2.6) and (2.7) over a doubly periodic domain fi. Let u'j 
{j — 1, 2) be the unique periodic solutions over Q, of the equations: 

J^u^^dx = 0, J -1,2, 

and set u'^ — — It is clear that, for C 7^ 0, the condition (2.10) ensures the 
smoothness of the function e"3 . 

Set Uj = v'j + Vj {j = 1,2). Then the system of equations (2.6) and (2.7) subject to 
periodic boundary conditions become: 

Avi = Ai (m\A\h"'<+"'''' + |5|2e«?+«§+^i+^2 + |c|2g«i+^i-^2 _ j + ^^pj^^ (4_2) 

Av2 = A2 (^|s|2e"?+"°+''i+^2 - |C|2e"3+''i-^2 _ + i^.^^, (4.3) 
^1,^2 are doubly periodic on dfl. 
The problem above defines the Euler-Lagrange equations of the functional: 
7(^1,^2) = 

2Ai 2A2 

- (Ci - ^) ^1 - (^2 - ^] V2] dx, V, e H{n), J = 1, 2. (4.4) 



Again, on the basis of the Moser-Trudingcr inequality (cf. [2, 15]), we check that / G 
C^{H{Q) X H(Q)). It is also easy to check the weak lower semi-continuity of I. 

On the other hand, integrating the equations (4.2) and (4.3) over Q, we obtain the 
natural constraints 

/■ (^|^|2e-S+-i + |ij|2e«?+«e+-i+'^.) = l{Ci + ^2m-27r(^ + ^] 
^ 2 / 2 V '^i / 

= 771 >0, (4.5) 

/ f-|A|2e-"?+-- + |C|V^+---) dx = l(^,-^,)\n\-2n(^-^] 
^ 2 / 2 V '^i / 

= 772 > 0. (4.6) 

Thus we deduce the necessity of the conditions (2.13) and (2.14) for the existence of a 
solution once we take into account of (2.9). 

Now assume that the conditions (2.13) and (2.14) are fulfilled. We show that (2.6) 
and (2.7) will have a solution. 
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To this purpose we use the decomposition (3.6) to obtain the following: 

Jn 

Jn 

> \B\'^\n\e^'+^^ -rii{v^ + V2) + \Cf\n\e^^-^^ -r]2{Vi-V2), (4-7) 

where again wc have used the Jensen inequality. 

At this point, by the non- degeneracy condition (2.9) and the fact that ?7i > and 
772 > (see (4.5) and (4.6)), as above we deduce that / is coercive and bounded from 
below in H{Q) x H{Q). By weak lower semi-continuity, we see that / attains its minimum 
value. Since / is strictly convex, the minimum point gives the unique critical point of /, 
and the first part of Theorem 2.2 is proved. 



5 Proof of Theorem 2.2: planar case 

In this section, we study the equations (2.6) and (2.7) over the entire plane subject to 
the topological boundary condition (2.12). As usual, || • ||p denotes the norm of L^(M^) and 
Bj. denotes the disk in centered at the origin with radius r > 0. We will continue to 
assume the condition (2.10) imposed on the sets of vortex points. Our goal is to estabhsh 
an existence and uniqueness theorem as in the doubly periodic situation. 
To this end, we introduce the background functions: 



/ I _ |2 \ 



AT,- 



(5.1) 



*W = ^"Ef^^Tk^^. ^ = 1.2- (5.2) 

We now set 

H,{x) = hi{x)h2{x), H2{x)^^, (5.3) 

11 11 
= T-c/i(a;) + T-c/2(a;), G2{x) = — ^i(x) - -^g2{x), (5.4) 

Ai A2 Ai A2 

which are all smooth functions. 

Setting Uj = + Vj (j = 1,2), we see that the equations (2.6) and (2.7) subject to 
the boundary condition (2.12) become the following boundary value problem over R^: 

Avi = mAi|A|2(/i™e™^i - 1) + Xi\B\\Hie'''+''^ - 1) + Xi\CWH2e'"-''^ - 1) + gi, 
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(5.5) 

Av2 = A2|S|2(i/ie^i+^^-l)-A2|C|^(i/2e''^-''2-l)+^2, (5.6) 
vi ^0, 0, \x\ oo. (5.7) 

As in the doubly periodic situation, we pursue a solution of the problem (5.5)-(5.7) 
by a minimization procedure. For this purpose, we observe that the problem enjoys a 
variational formulation with the action functional: 

I{vi,V2) = 

T^W'^Ml + ^\\Vv2\\l + \A\' [ {hT[e^'^-l-mv,] + [hT-l]mv,)dx 

IA\ IA2 Jm? 

+\B\^ f (i/i[e"i+"^ - 1 - (^;i + V2)] + [H, - l][v^ + V2]) dx 
+ |C|2 / {H2[e^'-''' -l-{vi- V2)] + [H2 - 1]K - V2]) dx 
+ / ( T-S'i^i + ^5'2^2 J dx, (5.8) 



Ai"^ A 



2 



where Vi,V2 G W^''^ 
Notice that. 



0<hj<l, 0<Hj<l; {1- h]'),{l- Hj) e LP{R'^),Wp> 1; 
5,eLi(R2)nL~(R2); J = 1,2. 



(5.9) 



So the functional (5.8) is well defined and of class in W^^^{R'^) x W'^'^iR"^). 

Moreover, exactly as for the Abelian case, elliptic estimates and standard regularity 
theory ensures that a critical point of / in Ty^'^(M^) x W^'^{R^) is a solution of (5.5)-(5.7). 
Furthermore, it can easily be checked that / is weakly lower semi-continuous and strictly 
convex. So it allows at most one critical point. 

Fix A > max{Ai, A2} > and rewrite (5.8) as follows: 

Iiv-i,V2) = 

+\A\^ [ (/i^^ie"*"! - 1 - mvi] + [/i^ - l]mvi) dx 

iR2 

+ \B\^ f {H^[e'''+^^ -l-{v^ + V2)] + [H^ - l]bi + V2]) dx 

Jr2 

+ |C|2 / {H2[e^'-''' -l-{v^- V2)] + [H2 - l][vi - V2]) dx 
+ / ( 7- S'l^'i + T- 5-2^^2 ) dx 

Jr2 \Ai A2 / 
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+ / ( -r5'i«'i + -^92V2 ] dx, (5.10) 



Ai ' " A 



2 



where 

^o(^) = ^11^^112 + "^' (x^-x) Vrj^^(^r[e^-l-^] + K-lHdx, (5.11) 
^i(^) = 1^11 + 2X13]^ [ {H,[e''-l-v] + [H,-l]v)dx, (5.12) 
Mv) = ^l|V^^||2 + 2A|Cp / (iJ2[e"-l-^;] + [i/2-l]t')dx. (5.13) 

These functionals may collectively be expressed in the form: 

J{v) = ^||V^;||^ + [ {HoAe" - 1 - v) + (i/o,^ - 1)^;) dx, (5.14) 

where a > is a constant and 

«o,.(x) = J]ln(^-^^^j, i/o,. = e-'^ (5.15) 

with Pi, ■ ■ ■ ,Pn G fixed points repeated accounting to their multiplicities. We have 
< Ho,^ < 1 and 1 - Ho,^ G L^{M.'^). Besides, for 

i?o >3max{|ps||s = l,---,A^}, (5.16) 

there exists some cq = Co(-Ro) > 1, such that for any /x > 1 we have: 

/^o,m(^) > r^, \x\>Ro, (5.17) 



Q<Ho,^{x) < ^, \x\<Ro, (5.18) 



/ HoAx)dx > (5.19) 
Jbr^ Co// 

Lemma 5.1. For the functional defined in (5.14.), there is a constant > 1 such that 
for II > /la there holds: 

J{v) > |||VH|^ + f Ik-ll + ^\M\' -C„ ve W'^'{R% (5.20) 

where = max{-;;,0} and v~ = max{— f,0} are the positive and negative parts of v, 
respectively, bo > is a constant independent of /j, and > is some constant depending 
on 11. 
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Proof. Since J{v) = J{v~^) + J{v ), we can split the proof into two cases 
Case 1. V <0 a.e. in R^. 
Recall the Sobolev inequality 

[ v^dx<2( [ v^dx]([ \Vv\'^dx], veW^^^R^). 

Thus, we have 

2 / /• I I \ 2 



/ n ^ I h2 / (^^ + '"^) 



< 2 

Inserting (5.21) into (5.22), we have 



Jr2 U + \^\) 
Besides, from the well-known inequality: 



1 - e-* > — ^, > 0, 



we deduce that. 



e-* 



l+t = /'(l-e-)d.>/^d.>^. W>0. 



Hence, for v < 0, we obtain the lower bound: 

> e^-l-w 

> 1 



21 + |'f;| 

Consequently, using (5.23), we have 



^ ^ - 2" "2 2 Ul + \v\ 



1,, ,,o a /" V 



,2 



- 2" 2 7m2 (1 + 1^^1)' 

> 7l|V^ll2 + 7 l|V^ll2 + 
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Therefore, we can minimize the function 



Vit) = t+-^^ t>0, (5.28) 



to find ip{t) > 2a\\v\\2 — 1, and conclude: 



1 rv 

> ^||V^;||^ + -||^;||2-1, ^;<0. (5.29) 



Case 2: v>0 a.e. in R^. 
First, note that 



/ (i^o,^(e^ -l-v) + {Ho,^ - l)v) dx> [ (iHoy + {Ho,^ - l)v] dx, (5.30) 
and from (5.17), we get 

/ (lHo,y + {Ho,, - l)v) dx 

J\x\>Ro \^ J 



\x\>Ro 

- 2^^"^"^'('^'\-^«o) ~ 11-^ ~ ^0,^IU2(M2\BHo)|k||L2(ffi2\Bfio) 



for a suitable constant Ci^^ > 0. 
On the other hand, if we set 



M = [ vdx, vi^v-M, (5.32) 



we find: 



^ {l^oy + {Ho,, - l)v^ dx 

= 2 / ^0,1^^1 + / Ho,^vi dx + — Hq,^ dx 

/ Ho,,vidx+ / {Ho,^-l)dx] 

If f f 

> 2 / Ho,f,vl dx+ Hq^^vi dx + — i^o.^ 

- I / HQ,^dx\ I I {Ho,,Vi + Ho,^ - 1) dx 
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If of f 

> - / Ho,,j,v^ dx+ Hq^ixVi dx + — Hq^^ dx 



/ Hq^iu dx 



Hq,i^Vi dx 



Ro 



Ho,ij,Vi dx 



Ro 



j (l-Ho,^)dx + f (l-Ho,^)dx 
[ Ho,^vidx] \ 2 f 



M2 

> — / i/o,/.dx + 



dx 



-1 



> 



M2 

T 

1 

'2 



Ho^fj^dx - c^\\vi\\l'2^Br^) 



Ro 



Ho,f, dx 



Ro 



Ro / 



for some suitable constants c^, C2,ij, > 0, after using the inequality 



Ho^^Vi dx> \ i^o.M da; 



Ho 



Rq 



Ho,u.vi dx 



Ro 



As a consequence, using (5.31), for v > 0, we have 



Ro 



-a 



+ 



a 



for some constant Cs,^ > 0. Now notice that 

-1 



Ho,^l dx 



Ro 



Hl^,dx < ||i/o,Ml|L-(Bflo) 



< 



Co 



Ro 



AT' 



in view of (5.18). On the other hand, by the Poincare inequality, we have 

lki|li2(B,^)<aoi?^||V^;i||i.(5^^), 
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which leads us to obtain: 



(5.38) 



Thus, by (5.18), (5.19), (5.35), (5.37), and (5.38), we can find constants bo,bi > 0, 
depending on Rq but independent of //, such that 



A* 



(5.39) 



Summarizing Cases 1 and 2 above, i.e., (5.29) and (5.39), we can find /^a > 1 sufficiently 
large such that (5.20) holds whenever fx > fia- 

We are now ready to prove that the functional / defined in (5.8) or (5.10) is coercive. 
In fact, applying Lemma 5.1 to the functional J^, £ = 0, 1, 2, defined in (5.11)-(5.13), we 
can find constants a, 6 > (independent of fi) and > 1, such that for > sufficiently 
large and N — max{A^i + N2, mNi}, we have 

I{vi,V2) > a {WVviWl + ||Vt;2||^ + |A|||t;i_||2 + \B\\\{vi + V2)-h + \C\\\{vi - V2)-\\2) 

+^ + miv, + V2h\\i + \cn{v^ - V2h\\i) 

A* 



+ / ( J^9ivi + J^92V2 ] da; - c^. 



M2 



(5.40) 



Recalling (5.4), for 5 7^ 0, C 7^ 0, we can estimate 

/ ( T~S'if^i + T-9'2f^2 ) dx = J / Gi{vi + V2) dx + ^ G2(^^l -^^2)dx 

1 „ 1 

> -^ll^^llhlPl +'"2112 - 211^21111^1 - ^2112 

Q 

> (Ikl + 1^2112 + ll^^l - ^^2||2) 

A* 

> -- (||(^;i + ^;2)-||2 + 11(^^1 - ^^2)-||2) 

A* 

-4r (11(^1 + ^2)+y + \\{vi - V2U\\l) - C,,„(5.41) 
A* 

where £ > is arbitrarily small. Incidentally, also notice that in the case C — but 
A 7^ 0, S 7^ 0, we can estimate 



^ givi + -^g2V2] dx ^ [ 



Ai 



A\ A2 



vi + ^52^1 + '"2]^ dx 



> — (llt'ilh + 11^1 + ^^2112) , 
A* 



(5.42) 
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and argue as before to arrive at a lower bound like that in (5.41). Similarly, in the case 
B = but ^4 7^ 0, C 7^ 0, we can estimate 



1 1 » ^ 

■T-fi'i^^i + ^92V2 J da; 

Ai A2 



1 1 
Ai A2 



^^1 - J^92[vi - V2] ) dx 



> --(||^^i||2+ Iki -^^2112), (5.43) 
/i 



so that a similar conclusion holds. In view of these results, we see that (5.40) allows us to 
obtain some constants a, b, c, Co > such that the following coercive lower bound holds, 

I{vi,V2) > a(||Vvi||^ + ||Vt;2||^) +6(||^;i_||2 + 1^2-112) 

+c(lki+||2 + 11^2+11^) -Co, vi,V2eW'^^{m'). (5.44) 

Therefore, by virtue of (5.44), we can use standard arguments to obtain the existence 
of a global minimum point for I as its unique critical point in W^''^{R^) x 1/F^'^(]R^), 
provided at least two of the constants A, B, and C do not vanish. Hence, the following 
slight extension of part (ii) of Theorem 2.2 holds: 

Theorem 5.2. Consider the equations (2.6) and (2.7) over the full plane such that the 

vacuum constraint (2.8) is satisfied with at least two of the constants A.B, and C not 
vanishing. Then there always exists a unique solution (^1,^2) satisfying the boundary 
condition Ui,U2 as \x\ 00. 

6 Proof of pointwise estimates 

In this section, we analyze equations (2.6) and (2.7) under the vacuum constraint (2.8) 
so that they take the forms: 

Aui = Ai (m|A|2[e"^"i - 1] + |sne"^+"^ - 1] + iCHe"^""^ - 1]) +47r^5p,_^, (6.1) 

N2 

Au2 = A2(|5|2[e"^+«2-l]-|Cne"i-"^-l])+47r5^(5p,,,. (6.2) 

Actually, after some simple manipulation, it is easy to rewrite (6.1)-(6.2) as follows: 
Aui = Ai (ml^l'ie'""! - 1] + [\B\^ + \C\^][e''' - 1] + [\C\^ - \B\^]e'''-'''[l - e"^] 

Ni 

+2|5|V^[coshM2-l]) +47rJ]5p,_^, (6.3) 

s=l 

Au2 = A2 (2|S| sinhii2 + [|C|2 - |S|2]e"i-"^[e"^ - 1] + [\Cf - \Bf][l - e"^]) 

N2 

+47rJ]5p,,.. (6.4) 
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Prom (6.3)-(6.4), by the maximum principle and with the aid of a continuity method, 
it is possible to show that for \B\ < \C\ every periodic solution (1*1,^2) of the equations 
satisfies mi < and M2 < 0. 

Here we focus on the physically relevant situation where 



\B\ > \C\ i.e. ,^2 > 0. 



(6.5) 



In the following, the solution (mi,M2) considered will either be periodic or satisfies (6.3)- 
(6.4) over together with the boundary condition (2.12). 

Also notice that in the periodic situation considered over the cell domain Vt and under 
the vacuum constraint (2.8), the necessary and sufficient conditions for existence take the 
form: 



m\A\^ + 2\B\^ > 
mU|2 + 2|CP > 



47r 




N2 


W\ 


Ui 








W\ 


Ui 


A2 



(6.6) 
(6.7) 



and provide explicit upper bounds on the vortex numbers A^i and N2 in terms of the 
size of fl. In addition, notice that (6.6)-(6.7) allow Ni — N2 — (i.e., the absence of 
vortices) which yields only the trivial solution ui — U2 — 0, consistently with the physical 
interpretation. 

So from now on we assume A^^i > and start by showing the following: 



Lemma 6.1. Let X2 > Xi, \B\ > \C\, and assume that 
Then 



m2UP<2(^ 



CP. 



Ui U2 

— ± — < 0, ana in particular Ui < 0. 

Ai A2 



(6.8) 
(6.9) 



Remark 6.2. It is important to notice that our assumption (2.15) compares to (6.8) as 
follows. If < m < 2, then (2.15) is stronger than (6.8) and implies it. If m > 2, then 
(2.15) and (6.8) are one and the same. 



Proof. Observe that, if by contradiction we assume that there exists Xq: 

—ui{xo) - —U2{xo) = max —ui - —U2 > 0, 
Ai A2 \Ai A2 / 

then xo ^ • • ■ ,Pi,nA ^ {P2,i, ■ ■ ■ ,P2,N2} and it satisfies 

> a(^-^] {xo) = m\A\'' (e"^"i(^o) - l) + 2|C|' (^eMxo)-u2{xo) _ 
\Ai A2 / 

12 2\C 



(6.10) 



> Ai(m2|>l|2 + 2|C|2) 



Ai A2 / A2 



U2{Xo). 
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Thus, from (6.8), we get U2{xq) > 0. As a consequence, ui{xo) > {Xi/X2)u2{xq) > 0. 
Therefore, there exists Xi ^ {pi,i, ■ • • ,Pi,Ni} such that 

ui{xi) — maxrti > 0. 

But then from (6.3) we readily see that U2{xi) < and we conclude with 

ui{xi) ^ ^i(a:i) _ ^2(3^1) ^ uijxo) _ Uijxo) ^ Ui{xo) ^ Ui{xi) 



Ai Ai A2 Ai A2 Ai 

That is, ^2(3^1) = U2{xo) — and by virtue of (6.3) again we have 



Ai 



/ N n «l(^o) U2{Xo) 

msixui — Ui[Xi) — [j — — : : = max 



Ai 



Consequently, 



A(^-^) > {m'\A\'-2\C\' 



A2 



Ui U2 

Ai A2 



> 2|C|%(^-^). 



Ml _ U2 

Ai A2 

So, we can use the strong maximum principle to conclude that 

Ui _ U2 

Ai A2 ' 

which is certainly never allowed by (6.1)-(6.2), unless ui = U2 = 0. That is, Ni = N2 = 0. 
Hence for iVi > 0, there holds 

(6.11) 



Ai A2 



On the other hand, from (6.11) we also get that, necessarily, Ui < 0. Indeed, if, 
as before, we assume that there exists Xi ^ {pi,s | s = l,---,A^i} such that Ui{xi) — 
max Ml > 0, then U2{xi) < and wc get a contradiction as follows: 



Q ^ UljXi) ^ UijXi) _ U2{Xi) ^ ^ 



Ai 



Ai 



A2 



In turn, we see that 



Ul ^U2 

Ai A2 



N2 



m|A|'(e'""^ -l) + 2|S|"(e 



-11 



\^ f 47r v-^ . 



> {m^\A\^-2\B 



Ai A2 



=1 



Ui U2 



s=l 



> 2|isrA.|^ + 5l. 
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where the last inequahty follows by the observation that 



m'\A\'-2\B\'iy-l] <m'\A\'-2\C\'(^-l] <0 and m<0. 



Whence, 



Ai A2 



as it follows easily by the maximum principle. 

Next we show how the stronger assumption (2.15) (when < m < 2) actually implies 

Ui±U2 < 0. 



To this purpose, observe that 

A{ui + U2) = Ai ( m|^|2[e"^"i - 1] - 



Ai 



- 1 



\C\^[e"' - 1] 

Ni N2 



+(Ai + A2)|Sr(e-+- - 1) +47r J^^,,, + Y.^p.,s , (6-12) 



A(iii - U2) = Ai m|^|2[e"^"i - 1] 



,s=l s=l 
|^|2jg«l+«2 _ 1] 



+(Ai + A2) |Cr(e"-"^ - 1) + 47r 5,,, - ^ 5,,,, 

\s=l s=l 

Therefore, if by contradiction we assume that there exists 

^ ^ ■ ■ ■ ,Pi,Afi} : (wi + «2)(^) = max(Mi + U2) > 0, 
then ui{x) > —U2{x), and we find: 

'A2 



(6.13) 



> A{ui + U2){x) > Ai m|^ne™"i(^) - 1] - 



Ai 



- 1 



|C|2[e2"i(^) - 1] 



In case < m < 2, then realizing that ui < (by Lemma 6.1) and by (2.15), we obtain 
the desired contradiction as follows: 



m\A\'^{e^M^)-l)- 



h 

Ai 







icr) 









|C|2^ (e2"i(*)-l) > 0. 



In case m > 2, then (6.8) and (2.15) coincide and in this case, we arrive at the desired 
contradiction, by observing that under (6.8), the function 



fit) = m\A\\e"^' _ 1) _ ( ^ _ 1 ) \C\\e'' - 1) 
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satisfies /(t) > 0, Vt < 0. 

Similarly, we show that Ui — U2 < 0. Since if by contradiction we assume there exists 

^ ^ {Pi,i, • ■ ■ ,Pi,iVi} \ {P2,i, • • • ,P2,iV2} : (^^1 - U2){x) = max(iii - U2) > 0, 
then ui{x) > U2{x) and we get 

> A(mi - U2){x) > Ai (^m|A|2[e'""i(*) - 1] - ^ - 1 |S|2[e2"i(^) - 1]^ . 

At this point, recalling that \B\ > \C\, we reach a contradiction arguing exactly as above. 
To conclude the proof of Theorem 2.2, it remains to estabhsh (2.17) and (2.18). In the 

periodic case, they just express (4.5) and (4.6) once we introduce the vacuum constraint 
(2.8) and use (2.16). In the planar case, they follow exactly as in the Abclian case (e.g., 
see the proof of Proposition 3.2.4 in [41]). Indeed, multiplying (5.5) by the test function 
Xr{x) = x(D with X £ C^{B2) satisfying % = 1 in 5i and integrate over R^, to obtain: 

/ {m|^|'(l - e"*"^) + W{1 - e"^+"^) + |C|'(1 - e"^""^)} Xiidx 

^ "T ( / 9iXR<ix-[ AviXRdx 



=-(/ 



smce 



and 



giXR'^x- I ViAxijdx 1 ^> 47r^ asi^^^oo, 

M2 / '^l 



lim / giXR'^x= I g^dx = Av:Ni, 
Jm2 Jr2 



viXr dx 



iR<\x\<2R)) (^J |Ax| dx^ ^0 asR^ 00. 



Thus, we deduce that (1 - e"*"!), (1 - e"i±"2j ^ L\R'^) and 

4 

l-«2 II _ _ 



m|A|2||l -e"™^||i + l^nil -e"i+"^||i + \C\^\\1 - e'^'-^'^Wi = (6.14) 



Using (5.6), by a similar argument, we get 

|SH|l-e"i+"2||i-|CH|l-e"i-"^||i = ^ / g2dx^^^. (6.15) 

A2 Jk2 A2 

So (2.17) and (2.18) are established. 

In concluding this section, we mention that, by standard methods, it is not hard to 
show that, under the given assumption, ui ± U2 approach zero at infinity exponentially 
fast. 
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7 Degenerate cases: doubly periodic solutions 

In this section, we focus on the problems over a doubly periodic domain fi. 

(i) A 0, B — C — 0. In this case the system decouples and we reduce to analyze the 
single equation: 

Ni 

Am = X^{m\A\\^--i,)+4nJ2^.M, (7-1) 

s=l 

which is the well-known Abelian vortex equation, that admits a unique solution if and 
only if 

47r-i<6|^^|- (7.2) 

(ii) A = 0, B ^ 0,C = 0. The equations (4.2) and (4.3) now become 

A^;i = Ai (^|s|2e"S+«§+^i+^2 _ + ^^^^^ (7 3) 

Av2 = A2(|S|V?+"^+''i+''^-6) + pAr2, (7.4) 

and no longer we need to assume (2.10). Since (4.6) gives us 

772 = 0, (7.5) 

then 

which implies that V2 — Xi^X2Vi up to a constant. The equations are turned into a single 
one, and a unique solution exists if and only if (7.2) holds. 

(iii) A — 0, B — 0,C 0. The governing equations now read 



Av^ = Ai(|C|V3+''i-''^-6)+^A^i, (7.7) 
Av2 = -A2 (|C|V3+"i-"2+e2) +^iV2. (7.8) 



Since (4.5) implies 771 = or 

Ai A2 



(6 + 6)l^^|-47r(^ + ^) =0, 



A(^ + ^)=0 (7.9) 



we have 

VAi A2 

and, thus, V2 = — Aj^^A2fi up to a constant, and we arrive at the same conclusion as in 
the case (ii). 
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(iv) A 7^ 0, S 7^ 0, C = 0. Again, we no longer need to require (2.10), as the equations 

are 

A^;i = Ai (m\A\^e"'<+"'''' + |5| V?+"°+^i+^2 - ^i) + ^A^i, (7.10) 



Av2 = A2(|S|V?+"^+-^+-^-6)+p7V2, (7.11) 
with associated action functional 

Jn 1 2Ai 2A2 J 

-ViiUi + V2) - miv^ - V2). (7.12) 

On the other hand, with C = 0, the constraints (4.5) and (4.6) are refined into 

/"!^|A|2e"^"?+"»^idx = 772 >0, (7.13) 
Jn 2 

f |^|2gn?+u0+^i+.2^^ ^ rji-ri2 = V3>0. (7.14) 
Jn 

Thus, we can substitute 

Vi{vi + V2) + V2{vi - = 27^21^1 + ^3(1^1 + V2) (7-15) 

into (7.12) and follow the proof of Theorem 2.2 to show that the functional (7.12) has a 
unique critical point in H[Q) x H[il). 

(v) A 7^ 0, S = 0, C 7^ 0. The equations are 

Avi = Ai (^m|A|'e"^"?+"*^^ + |C| Va+^'i-^^ - + ^AT^, (7.I6) 

Av2 = -A2(|C|V3+''^-^^+e2)+piV2, (7.17) 

with associated action functional 

I(v,,V2) = / (t^IV^^iT + 7^|VT;2p + |Ape"*«?+--^ + \C\\-°s+-^-A dx 
J a L2Ai 2A2 J 

-Vi {vi + ^^2) - ^2 (vi - V2) • (7-18) 
When B — 0, the constraints (4.5) and (4.6) become 

[ ^\A\^e"'<+'^'" dx = 771 >0, (7.19) 
Jn 2 

|^|2g«o+t,i-t,2 _ ^2 - r^i = 7^4 > 0. (7.20) 
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Therefore, we may insert 



ViiHi + II2) + V2{vi - V2) = ^ViV.! + ^4fe - V2) (7-21) 

into (7.18) and adapt the proof of Theorem 2.2 to estabhsh the existence and uniqueness 
of a critical point of the functional (7.18). 
In summary, we can state 

Theorem 7.1. In the degenerate cases when B or C vanishes (we need to require (2.10) 
only when C ^ 0), the existence and uniqueness of a solution to the master equations 
(2.6) and ( 2. 1) hold under the following necessary and sufficient conditions. 

If A 0, B — C — 0, the equations are reduced to (7.1) and the condition is (7.2). If 
A = 0, S 7^ 0, C = or A = 0, S = 0, C 7^ 0; the equations are also reduced to (7.1) and 
the condition reads 



along with 



or 



respectively. 

IfA^O,B^O,C — 0, then the condition reads 



/iVi No\ No 
47r - ^ j < (a - 6)1^^1, 47r^ < (7.25) 

IfA^O,B = 0,C^O, then the condition states 

/A^i iVoX No 
47r(^ + ^J < (6 + 6)1^^1, 6|^^|<47r^. (7.26) 

In all cases, the solutions may be constructed by a direct minimization method. 

It is interesting to note that the pair of conditions (7.23) and (7.24), and (7.25) and 
(7.26), are interchangeable under the correspondence (iV2,^2) —(-^2,^2)- 



8 Degenerate cases: planar solutions 

We now consider the degenerate cases concerning the planar solutions of the equations 
(2.6) and (2.7) subject to the boundary conditions (2.12). We have already shown exis- 
tence and uniqueness when at least two of the constants A, B, and C do not vanish. Here 
we consider the other degenerate cases. 
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(i) A ^ 0,B = C = 0. In this case, we do not need (2.10), and the condition (2.8) 
gives .^2 = and the system decouples into an equation for Ui as that in the Abehan Higgs 
model (1.6) and a linear equation for U2, given as follows: 

Au, = X,{m\A\'e"^^^-i^)+47TJ2Sp.M, (8-1) 

N2 

Au2 = A7rY,5,,^^{x). (8.2) 

s=l 

Obviously, the natural boundary condition for (8.1), (8.2) only involves Ui and requires 
Ui{x) as I a; I oo. Thus the existence and uniqueness of a solution ui < follows 
by the same analysis of the Abelian Higgs model (1.6). While a solution of (8.2) with 
logarithmic growth may be obtained via the fundamental solution. 

(ii) A^O,B y^O,C ^0. Prom (2.8), we see that = ^2 = C = Thus the system 
becomes 

Am = Ai(|5|Vi+'^^-e)+47r^5,,,^(x), (8.3) 

s=l 

N2 

Au2 = A2(|S|Vi+"^-e)+47r^(5^,,^(x). (8.4) 

Now the natural boundary condition involves only Ui + U2 as follows: {ui + U2){x) — )■ as 
— 7> 00. With the background functions defined in (5.1) and (5.2) and setting Uj = Uj+Vj 
{j = 1,2), we may recast the above equations into 

Avi = Ai (^|5|2e«?+«e+^i+^2 - ^) + (8.5) 
Av2 = A2 V?+"°+^^+^^ - + ^2. (8.6) 

Again, we no longer need to require (2.10). Moreover, (8.5) and (8.6) can be easily 
handled as the Abelian model (1.6) once we introduce the known w — vi + V2 that must 
satisfy: 

Aw = (Ai + X2)\B\\H,e^ -l)+gi + 92, .g 
ty — >■ as |x| — >■ 00, 

where the function is as given in (5.3). 

We know the existence and uniqueness of a solution w e l^-'^'^(R^) of (8.7) that also 
satisfies: i?ie"' < 1, ty — > exponentially fast at infinity, and 

|i;r/(l-g.e")d.= ^^f' + ^-' . (8.8) 
Jr2 Ai + A2 

At this point, we can use V1 + V2 = w m. (8.5)-(8.6) and linear elliptic theory to find Vi and 
V2 (unique up to an additive constant) and see that in general they admit a logarithmic 



24 



growth at infinity at a rate determined by the integral value of the right-hand side of (8.5) 
and (8.6) respectively. 



From (8.8), we easily compute: 

— (xi\B\^f {Hie"" -l)dx + 47rNi 
27r V 



2A1A2 fNi_N2 
Ai + A2 \ Ai A2 



as \x\ 00. (8.9) 



and obtain: 

-i(^) = lrTfe(t-t)inN + o(i), 
-2(^) = -|r^(t-t)i-N + o(i)' 

For details, see e.g. Lemma 1.1 and Corollary 1.1 of [7]. 

Thus, in this case in terms of the original variables Ui and U2, we can claim the 
existence of a solution for (8.3)-(8.4) satisfying: ui + U2 < 0, Ui + U2 ^ exponentially 
fast at infinity, and 



A1+A2 

Moreover, by the maximum principle, we see that: 





N2 
A2 


<0, 




N2 
A2 


>0, 


. Ni _ 
■ Ai 


N2 
Ai ' 


we c 



^)x 



W^^'^(]R^); and so both ui and U2 vanish at infinity and mi < 0, M2 < in R^. 

(iii) A — B — 0, C 0. This case requires (2.10) and the natural boundary condition: 

{ui — U2){x) — >■ as |a;| — >■ 00. 

Such a case can be treated as above with the help of the new variable w = Vi — V2, 
which must satisfy a problem as (8.7) only with Hi replaced by H2, by |Cp, and 
gi + 92 by gi - g2- 

So, by similar arguments, also in this case the existence of a solution pair (ui, ^2) can 
be established satisfying: i^i — < in R^, — ^2 — >■ at infinity exponentially fast, 
and 

,2 [ n .™ - 47r(iV, - N2) 



\CY / (1 - e"i-"2) dx = 

Jr2 Ai + A2 

Finally, we remark that in all parameter regimes the following "flux quantization" 



formulas 



/47r7Vi 
(m|A|2e"*"i + |S| Vi+"^ + |C| Vi-"2 - Ci) = (8.10) 

Vi+"2 - |C| V^-"^ - 6) = (8.11) 

A2 
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are valid for solutions of (2.6)-(2.7) over a doubly periodic domain, or the full plane. These 
expressions are weaker than but cover those stated in (2.17) and (2.18) which involve the 
L^-norms of various quantities. They are exactly what have been used in literature for 
the computation of the tension of the non-Abelian vortex tubes. 
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